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We report experimental results on the properties of Bose- 
Einstein condensates in ID optical lattices. By accelerating 
the lattice, we observed Bloch oscillations of the condensate 
in the lowest band, as well as Landau-Zener (L-Z) tunneling 
into higher bands when the lattice depth was reduced and/or 
the acceleration of the lattice was increased. The dependence 
of the L-Z tunneling rate on the condensate density was then 
related to mean-field effects modifying the effective potential 
acting on the condensate, yielding good agreement with re- 
cent theoretical work. We also present several methods for 
measuring the lattice depth and discuss the effects of the mi- 
cromotion in the TOP-trap on our experimental results. 

PACS number(s): 03.75.Fi,32.80.Pj 



I. INTRODUCTION 

In a very short time after their first observation, Bose- 
Einstein condensates (BECs) have advanced from being 
mere physical curiosities to the status of well-studied 
physical systems Q. A host of diverse and interesting 
phenomena such as collective modes, quantized vortices, 
and solitons, to name but a few, have been extensively 
investigated and are now well understood Going 
from the harmonic potentials used in the experiments 
mentioned above to optical lattices ]7j constitutes a nat- 
ural extension of the experimental efforts to periodic po- 
tentials and has opened up new avenues for research. So 
far, experiments using periodic potentials have focused 
mainly on Bragg scattering [j§|-[lCfl and, more recently, 
on phase properties, involving such intriguing concepts 
as number squeezing 1 1 and the Mott insulator transi- 
tion 12 Some interesting work has also been done on 
superfluid properties of BECs in optical lattices IblJlF 



In the present work we report on experiments with 
Bose-Einstein condensates adiabatically loaded into one- 
dimensional optical lattices [p~6|Jl~7| ] . In particular, we look 
at the dynamics of the BEC when the periodic potential 
provided by the optical lattice is accelerated, leading to 
Bloch oscillations and L-Z tunneling. We then proceed 
to use L-Z tunneling as a tool for measuring the effects of 
the mean-field interaction between the atoms in the con- 
densate. The modification of the L-Z tunneling rate in 
the presence of interactions can be interpreted in terms 
of an effective potential, and we obtain good qualitative 
agreement with a recent theoretical study using this ap- 
proach nM. 



This paper is organized as follows: After briefly in- 
troducing some essential ideas and terminology used in 
the theoretical treatment of cold atoms in optical lat- 
tices (Sec. ||), we describe our experimental apparatus 
in Sec. |ni| . After presenting in Sec. |^ the results of pre- 
liminary experiments on the calibration of the lattice and 
on the effects of the condensate micromotion, we turn to 
the subject of Bloch oscillations in accelerated lattices in 
Sec. |v|. The following section VI deals with L-Z tunnel- 
ing and l eads on t o a d iscussion of mean-field effects in 
Sec. |VIl[ . In Sec. VIII we present our conclusions and 
an outlook on further studies, followed by an Appendix 
in which we discuss various parameters relevant to the 
description of our system as an array of tunneling junc- 
tions. 



II. COLD ATOMS IN PERIODIC STRUCTURES: 
BASIC CONCEPTS 

The properties of cold atoms in conservative optical 
lattices (i.e. far-detuned, so spontaneous scattering is 
negligible) bear a strong resemblance to the behaviour of 
electrons in crystal lattices in condensed matter physics 
and have, therefore, enjoyed increasing interest since the 
early days of laser cooling. There are a number of excel- 
lent review papers on the subject [|l9|,|o), so in the follow- 
ing we shall only briefly review some basic concepts and 
establish conventions and notations for the remainder of 
this work. 

The physical system we are considering consists of a 
Bose-Einstein condensate in a periodic potential created 
by two interfering linearly polarized laser beams with 
parallel polarizations. The potential seen by the atoms 
stems from the ac-Stark shift created by an off-resonant 
interaction between the electric field of the laser and the 
atomic dipole. This results in an optical lattice potential 
of the form 



U{x) = Uq sm 2 (ttx / d) , 



(1) 



where d is the distance between neighbouring wells (lat- 
tice constant) and 



u = (2/3)fir(j// )(r/A), 



(2) 



is the depth of the potential |l9[|, where / is the inten- 
sity of one laser beam, Iq is the saturation intensity of the 
87 Rb resonance line, T is the decay rate of the first excited 
state, and A is the detuning of the lattice beams from 
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the atomic resonance. If the momentum spread of the 
atoms loaded into such a structure is small compared to 
the characteristic lattice momentum pb — 2Tm/d, then 
their thermal de Broglie wavelength will be large com- 
pared to the lattice spacing d and will, therefore, extend 
over many lattice sites. The condensates used in our 
experiment have coherence lengths comparable to their 
spatial extent of ~ 10 fim, which should be compared to 
typical lattice spacings in the region of 0.4 — 1.6 ^m. A 
description in terms of a coherent delocalized wavepacket 
within a periodic structure is then appropriate and leads 
us directly to the Bloch formalism first developed in con- 
densed matter physics. As will be explained in the fol- 
lowing section, we could vary the lattice spacing d by 
changing the angle between the lattice beams. In this 
article, the lattice spacing d always refers to the respec- 
tive geometries of the optical lattices used. 

In a lattice configuration in which the two laser beams 
with wavevector kh are counter-propagating, the usual 
choices of units are the recoil momentum p rec — Tik^ = 
mv rec and the recoil energy E rec — h 2 k\/2m. In the 
case of an angle-geometry, with the angle 9 between the 
lattice beams (see Fig. ||), it is more intuitive to base 
the natural units on the lattice spacing d = k ^ si ^gy 2 ) 
and the projection k = ir/d of the laser wavevector k^ 
onto the lattice direction. One can then define a Bloch 
momentum pb — 2hn/d = mvs, corresponding to the 
full extent of the first Brillouin zone or, alternatively, 
to the net momentum exchange in the lattice direction 
between the atoms and the two laser beams. Possible 
choices for the energy unit are either the Bloch energy 
denned as Eb = ft 2 (2tt) 2 / md 2 , or an 'effective' recoil en- 
ergy E rec (9) = Eb/S, where the parameter 9 indicates 
the dependence on the lattice geometry |H) . As an intu- 
itive choice for the natural units for the lattice depths is 
the geometry dependent recoil energy E rec (9), we shall 
quote the lattice depth Uq in units of this scaled recoil en- 
ergy; for simplicity of notation we shall write E reci where 
it is understood that this always refers to the respective 
lattice geometries. In Sec. IV on the calibration of the 



lattice depth, we also use the parameter s = Uo/E rec . 
Throughout the paper, velocities and momenta will be 
quoted in units of the Bloch momentum ps and Bloch 
velocity vb- 

In the tight-binding limit (/To ^ 10 i? rec ), the conden- 
sate in the lattice can be approximated by wavepackets 
localized at the individual lattice sites (Wannier states). 
This description is more intuitive than the Bloch picture 
in the case of experiments in which the condensate is re- 
leased from a (deep) optical lattice into which it had pre- 
viously been loaded adiabatically. In the present work, 
this description is only made use of in Sec. |^ and in the 
Appendix, where the Wannier states are approximated 
by Gaussian functions. 



III. EXPERIMENTAL SETUP 

Our apparatus for creating Bosc-Einstein condensates 
is described in detail in Briefly, we use a double- 

MOT setup in order to cool and capture 87 Rb atoms 
and transfer them into a time-orbiting potential (TOP) 
trap. Starting with a few times 10 7 atoms in the mag- 
netic trap, we evaporatively cool the atoms down to the 
critical temperature for condensation in ss 30 s, obtaining 
pure condensates containing up to 2 x 10 4 atoms. After 
condensation, the magnetic trap is adiabatically relaxed 
to mean trap frequencies Ptrap on the order of 20 — 60 Hz, 
resulting in a variation of the condensate peak density 
between 2 x 10 13 cm -3 and 10 14 cm -3 . 

The optical lattice was realized using two linearly 
polarized Gaussian beams (waist ~ 1.8 mm, maximum 
power 3 mW) independently controlled by two acousto- 
optic modulators (AOMs) and detuned by about 30 — 
60 GHz above or below the rubidium resonance line. The 
lattice constant d could be varied through the angle 9 
between the laser beams, as shown in Fig. [l]. Both hori- 
zontal and vertical optical lattices with various angles 9 
were realized. Furthermore, by introducing a frequency 
difference 8 between the two beams, the lattice could 
be moved at a constant velocity vi at — 2 s\n(e/2) ^ or ac ~ 



celcrated with an acceleration a = 



,16 



While 



2sin(fl/2) dt ' 

in the countere-propagating geometry lattice depths up 
to ps 2 E rec were realized, in the angle-geometry lattice 
depths up to ps 20 E rec could be realized. 
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FIG. 1. Schematic of the experimental procedure. The con- 
densate is loaded into the optical lattice, which can transfer 
momentum to it in units of pa = tuvb- The frequency dif- 
ference S between the lattice beams can be used to create a 
moving or uniformly accelerated lattice. 



IV. PRELIMINARY EXPERIMENTS 



In order better to understand the effect of the lattice 
on the condensate and to calibrate the theoretically cal- 
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culated lattice depth against experimental values, we per- 
formed a series of preliminary experiments in conditions 
in which we expected mean-field effects to be negligible, 
i.e. either with the condensate in a weak magnetic trap 
(after adiabatic relaxation), resulting in a low condensate 
density, or by loading the condensate into the lattice af- 
ter switching off the magnetic trap. In the latter case, 
for horizontal lattice configurations the condensate was 
in free fall after switching off the TOP-trap, limiting the 
interaction time with the lattice to 10— 15 ms. In all ex- 
periments, the condensate was observed after 10 — 20 ms 
of time-of-flight by flashing on a resonant imaging beam 
for 20 /j,s and observing the shadow cast by the conden- 
sate onto a CCD-camera. 
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FIG. 2. Rabi oscillations of the condensate in an optical 
lattice. Shown here is the fraction Ni/N of condensate atoms 
in the first diffraction order as a function of time. The lat- 
tice (in the counter-propagating configuration) was moving 
with a constant velocity v = \vb- From the Rabi period 
TRabi ~ 260 /is one finds a lattice depth Uo « 2E rec . 



As pointed out in p9| , an alternative way of interpret- 
ing the observed Rabi oscillations in this kind of exper- 
iment is to consider a two-photon Raman coupling be- 
tween the two momentum states, whose energies differ 
by AE — 2h 2 k 2 /m = AE rec . The coupling is resonantly 
enhanced if the frequency difference S between the two 
lattice beams matches AE/h, i.e. if the lattice velocity 
viat — \ S = ^ = \vb-, as before. The two-photon Rabi 
frequency for the beam intensities corresponding to a lat- 
tice depth Uq can be easily calculated and, again, gives 

n R = u /2h. 
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FIG. 3. Calibration of the lattice depth by measuring the 
side-peak populations. We calculated the lattice depth from 
the mean value of the plus and minus first order side-peak 
using Eqn. ^. By varying the intensity of the lattice beams, 
we performed measurements of s exp for various values of Scale- 



2. Analysis of the interference pattern 



A. Calibration of the lattice 

1. Rabi oscillations 

If we abruptly switch on an optical lattice moving at 
a speed \vb-, then in the band-structure picture the con- 
densate finds itself at the edge of the Brillouin zone where 
the first and second band intersect at zero lattice depth. 
Raising the lattice up to a final depth of Uq opens up 
a band-gap of width AE = Uq/2 in the shallow lattice 
limit, and hence the populations of the two bands ]2~i| ] 
accumulate a phase difference A<f> = which results 
in the two populations getting back into phase (modulo 
27r) after a time TRabi = fr ■ ^ n ^ ne res t frame of the lab- 
oratory, one observes Rabi oscillations || between the 
momentum classes \p = 0) and \p = ps) in the shape 
of varying populations of the corresponding diffraction 
peaks observed after a time-of-flight (see Fig. [2]). From 
the oscillation frequency £Ir we could then calculate the 
lattice depth Uo — 2MIr, yielding results that fell short 
by about 20 — 25 percent of the calculated value. We 
attribute the discrepancy between the experimental and 
theoretical values mainly to the uncertainty in our laser 
intensity measurements and to imperfections in the lat- 
tice beam cross-section and polarization. 



If the depth of a stationary lattice is increased on a 
timescale comparable to the inverse of the chemical po- 
tential of the condensate (in frequency units), then the 
condensate can adiabatically adapt to the presence of 
the periodic potential |25|] . When the lattice has reached 
its final depth, the system is in a steady state with the 
condensate distributed among the lattice wells (in the 
limit of a sufficiently deep lattice in order for individual 
lattice sites to have well- localized wavepackets) . If the 
lattice is now switched off suddenly, the individual (ap- 
proximately) Gaussian wavepackets at each lattice site 
will expand freely and interfere with one another. (In 
this case a tight-binding approximation is more intuitive 
than a Bloch wave approach.) The resulting spatial inter- 
ference pattern after a time-of-flight of t will be a series of 
regularly spaced peaks with spacing vst, corresponding 
to the various diffraction orders, with a Gaussian popula- 
tion envelope of width ss Tit /ma, where a is the width of 
the wavepackets at the individual lattice sites. In particu- 
lar, Pedri et al. have shown p(| that the relative popula- 
tions P±\ of the two symmetric plus and minus first order 
peaks with respect to the zeroth-order central peak are 
given by P±i = exp(—Aw 2 a 2 /d 2 ) ^7|. For deep lattice 
wells, a/d can be found by making a harmonic approx- 
imation to the sinusoidal lattice potential about a po- 
tential minimum, giving a Gaussian width a = d/irs 1 / 4 . 
For the relatively shallow lattice used in our experiments 
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(s < 20), however, this approximation is not very accu- 
rate. Instead, we used a variational ansatz for a Gaussian 
wavepacket in a sinusoidal potential. The resulting tran- 
scendental equation can be solved numerically to yield 
a (see Ref. |2(| and Appendix). Alternatively, we can 
find an analytical expression for the lattice depth as a 
function of the measured side-peak population, giving 

a 16 p- 1 / 4 ( 3 ) 

Sexp ~ [ln(P ±1 )] 2 ±x ■ [6) 

This expression can be used directly to calibrate the lat- 
tice based on a measurement of the side-peak popula- 
tions. Figure || shows the lattice depth s exp as inferred 
from the above equation by measuring the populations 
of the zeroth and plus/minus first order diffraction peaks 
plotted against the lattice depth s ca ic calculated from 
the beam intensity and detuning, taking into account 
the losses at the cell windows (w 8 percent). A straight 
line fit gives s exp = (0.76 ± 0.1)s ca ; c , consistent with the 
results obtained by measuring the Rabi oscillation fre- 
quency. 
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FIG. 4. Calibration of the lattice potential with d = 1.2/im 
by measuring the tunneling rate. Plotted here as a function of 
the calculated lattice depth s ca i c is the number of condensate 
atoms Ntrapped remaining in the lattice at time ti a t = 10.1 ms 
after switching off the magnetic trap. Fitting a theoretical 
curve (see text) to the data we find that s exp w 0.76 x s ca i c , 
consistent with the results of the other calibration techniques. 



3. Tunneling 

If the lattice beams are arranged such as to create a 
periodic potential along the vertical direction, for deep 
enough lattice potentials the condensate can be held 
against gravity for several hundred milliseconds. As the 
potential is reduced below a critical depth the condensate 
starts tunneling out of the lattice. If this critical depth 
is fairly small (< 12E rec ), then the tunneling rate can 
be calculated from L-Z theory. If the condensate mov- 
ing with the acceleration a crosses the band gap AE at 
the Brillouin zone edge fast enough, the probability r for 
undergoing a transition to the first excited band is given 
by §§ 

r = e~ a ' /a , (4) 



with the critical acceleration 

a c = -^(AE) 2 . (5) 
An 

If the atoms can tunnel into the second band and, there- 
fore, effectively into the continuum, they will no longer 
be trapped by the lattice. Starting from this assump- 
tion, for the condensate accelerated by gravity with 
a = 9.81ms -2 we can calculate N tra pped the number of 
atoms that remain trapped after a time t\ at to be 

N trapped A^im(l ^ 3 (6) 

where ti at is an integer multiple of the Bloch period 
tb = v b/ 9 (see next section). Figure || shows the results 
of an experiment in which the condensate was loaded 
into a vertical lattice with lattice constant d = 1.2 /xm, 
after which the magnetic trap was switched off and the 
number of trapped atoms was determined after a time 
tiat — 10.1 ms = 26 tb as a function of s ca i c - Fitting 
the above equation for N tra pped(t) to the data, we found 
that the actual lattice depth was around 75 percent of the 
value calculated from the lattice beam parameters. This 
value agrees with those found using Rabi oscillations and 
side-peak intensity. 

B. Effects of the micromotion 

The time-orbiting potential trap used in our experi- 
ments is an intrinsically dynamic trap which relies on the 
fast rotation of the bias field to create an averaged har- 
monic trapping potential. It has been shown, however, 
that the atoms in such a trap perform a small but non- 
negligible micromotion at the rotation frequency Qtop 
(2tt x 10 kHz in our set-up) of the bias field ^9|. Al- 
though the spatial amplitude of this fast motion is ex- 
tremely small (less than 100 nm for typical experimental 
parameters), its velocity can be as large as a few mil- 
limeters per second. As the Bloch velocities vb of the 
lattices used in our experiments lie between 3 mms -1 and 
11.8 mms -1 (depending on the angle 0), the micromotion 
velocity component along the lattice direction can be a 
significant fraction of the Bloch velocity. This has two 
consequences: 

(a) If the interaction time of the lattice with the conden- 
sate is short compared to 27r/r2xop = 100 /is, the initial 
velocity of the condensate is undetermined to within the 
micromotion velocity amplitude. 

(b) If, on the other hand, the interaction time with 
the lattice is long compared with 27r/r2xoPj the conden- 
sate oscillates to and fro in the Brillouin zone and can, 
if the micromotion velocity is large enough, reach the 
band edge and thus be Bragg-reflected. Another possi- 
ble mechanism is the parametric excitation at Qtop of 
transitions to higher bands. 

All of these effects are undesirable if one wants to con- 
duct well-controlled experiments. In our setup, the mi- 
cromotion takes place in the horizontal plane and was 
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thus important when we worked with horizontal lattices. 
In order to minimize the detrimental effects of the micro- 
motion, we employed two techniques: 

(a) For short interaction times, we synchronized the in- 
stant at which the lattice was switched on with a given 
phase of the rotating bias field. This allowed us to en- 
sure that the lattice was always switched on when the 
condensate velocity along the lattice direction was ap- 
proximately zero. A small residual jitter, however, was 
still given by the sloshing (dipole oscillation) of the con- 
densate in the magnetic trap, which was especially crit- 
ical when we worked with large lattice constants (and 
hence small Bloch velocities). 

(b) For long interaction times, we phase-modulated one 
of the lattice beams synchronously with the rotating bias 
field and with a modulation depth that resulted in the 
lattice 'shaking' with the same velocity amplitude as the 
micromotion. In this way, in the rest-frame of the lat- 
tice the condensate was stationary (again, save a possible 
sloshing motion). 




FIG. 5. Effects of the micromotion in a horizontal 
counter-propagating lattice configuration, (a) shows a con- 
densate (after a time-of-flight of 15 ms) released from the 
magnetic trap with Vtrap = 25 Hz without the optical lat- 
tice. In (b), the optical lattice((7o ~ 2E rec ) was switched on 
for 7.5 ms with the magnetic trap still on. One clearly sees 
that the condensate is broadened along the lattice direction. 
In (c), we compensated the micromotion by phase- modulating 
one of the lattice beams (see text). The difference in intensity 
of the dark spot in (a) and (c) is due to shot-to-shot fluctua- 
tions of the number of atoms in the condensate. Also, in (c) 
a faint spot to the left of the central condensate can be seen; 
this corresponds to a part of the condensate having undergone 
Bragg reflection due to an initial sloshing. 

Figure || shows the effect of method (b) in which a 
condensate was loaded into a lattice (Uo ~ 2 E rec ) with 
a 1 ms ramp and then left in the lattice for 7.5 ms before 
the latter was suddenly switched off. When the micro- 



motion was not compensated by phase-modulating one 
of the lattice beams, the condensate appeared 'smeared 
out' when observed after a time-of-flight. Using the com- 
pensation technique eliminated this effect (Fig. || (c)). 




FIG. 6. Acceleration of a Bose-Einstein condensate in the 
counter-propagating geometry [d = 0.39 /im). In (a)-(f) the 
lattice parameters were Uo = 2.3 E rec and a = 9.81ms -2 , 
and the condensate was accelerated for 0.1,0.6,1.1,2.1,3.0 
and 3.9 ms, respectively. In (g), the condensate was accel- 
erated for 2.5 ms with the same lattice depth as above, but 
with a = 25ms -2 . In this case, a fraction of the conden- 
sate underwent L-Z tunneling out of the lowest band each 
time a Bragg-resonance was crossed. Note that the separa- 
tions between the spots vary because detection occurred after 
different times-of-flight. 



V. BLOCH OSCILLATIONS 

A. Theoretical considerations 

One of the most intriguing manifestations of the quan- 
tum dynamics of particles in a periodic potential are 
Bloch oscillations. Their theoretical explanation is based 
on the evolution in the band-structure picture of a col- 
lection of particles occupying a small fraction of the Bril- 
louin zone when the potential is switched on (meaning 
that in real space their wavefunctions extend over many 
lattice sites, which translates into temperatures well be- 
low the recoil temperature T rec = E rec /kB, see Sec. 0). 
If the lattice is switched on adiabatically, then all the 
atoms will end up in the lowest band. Accelerating the 
atoms by applying a force (real or inertial) to them will 
result in their being moved through the Brillouin zone 
until they reach the band edge. Owing to the effects 
of the periodic potential, at this point there is a gap 
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between the first and second band, and unless the ac- 
celeration is large enough for the atoms to undergo a 
L-Z transition (see Sec. VI), they will remain in the first 
band and thus be Bragg-reflected back to the opposite 
end of the Brillouin zone. In the rest-frame of the lattice 
this corresponds to the atoms' velocity oscillating to and 
fro between +^<xvb and —\olvb (where < a < 1, de- 
pending on the lattice depth) , whereas in the laboratory 
frame these Bloch oscillations manifest themselves as an 
undulating (or, for shallow lattices, almost stepwise) in- 
crease in velocity rather than the linear increase expected 
in the classical picture in which the atoms are 'dragged 
along' by the potential. The instantaneous velocity of 
the atomic wavepacket in the lattice frame can be calcu- 
lated from the slope of the first band at the corresponding 

1 dE{q) 
h dq 



quasi-momentum q, giving v 




1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 
FIG. 7. Bloch oscillations of the condensate mean ve- 
locity v m in an optical lattice. (a) Acceleration in the 
counter-propagating lattice with d = 0.39 /im, Uo ~ 2.3 E rec 
and a = 9.81ms -2 . Solid line: theory, (b) Bloch oscillations 
in the rest frame of the lattice, along with the theoretical pre- 
diction (solid line) derived from the shape of the lowest Bloch 
band. 



B. Experimental results 

The condensate was loaded into the (horizontal) opti- 
cal lattice with lattice constant d — 0.39 /im immediately 
after switching off the magnetic trap. The switch-on 
was done adiabatically with respect to the lattice vibra- 
tion frequencies ui v n, — 2E rec ^/s/h (valid for s ^> 1) 
by ramping up the lattice beam intensity over a time 
tramp ~ 100 [as. Thereafter, the lattice was accelerated 
with a = 9.81 ms" 1 by ramping the frequency difference 
5 between the beams. After a time t acce i, the lattice was 
switched off and the condensate was observed after an ad- 
ditional time-of-flight of 13 — 18 ms. Figure ^ shows the 
results of these measurements in the laboratory frame. 
The Bloch oscillations are more evident, however, if one 
calculates the mean velocity v m as the weighted sum over 
the momentum components after the interaction with the 



accelerated lattice, as shown in Fig. |7[ When the instan- 
taneous lattice velocity vi a t is subtracted from v m , one 
clearly sees the oscillatory behaviour of v m — vi at . This 
result is analogous to the first observation of Bloch os- 
cillations in cold atoms at sub-recoil temperatures [ jl9| . 
The added feature in our experiment is that by using 
a Bose-Einstein condensate released from a weak mag- 
netic trap, the spatial extent of the atomic cloud is suffi- 
ciently small so that after a relatively short time-of-flight 
(w 10 — 20 ms) the separation between the individual mo- 
mentum classes is already much larger than the size of 
the condensate due to its expansion and can, therefore, 
be easily resolved. The mean velocity is then calculated 
simply by counting the number of atoms in each of these 
classes (the dark dots visible in Fig. 0). 




FIG. 8. L-Z tunneling of a condensate in an optical lattice, 
(a) and (b): Mean velocity of the condensate after accelera- 
tion of the lattice to vb and 3vb, respectively, as a function of 
acceleration, (c) Mean velocity of the condensate after accel- 
eration of the lattice to 4.4vb as a function of lattice depth. 
In (a) and (b), the lattice depth was fixed at Uq — 2 E rec , and 
in (c) the acceleration a — 9.81 ms -1 . In (c), agreement with 
theory is expected to be somewhat less good because the fi- 
nal velocity of the lattice is not an integer multiple of vb (see 
text). 



VI. L-Z TUNNELING 

We investigated the Uq dependence for the L-Z tunnel- 
ing of the condensate into the second band when crossing 
the edge of the Brillouin zone, and therefore, effectively 
to the continuum, as the gaps between higher bands were 
negligible for the shallow potentials used in our experi- 
ments. As in Sec. [v|, we loaded the condensate into the 



G 



optical lattice after switching off the magnetic trap. The 
trap had been adiabatically expanded to a mean trap 
frequency v tra p ~ 20 Hz prior to switch-off, thus ensur- 
ing that the condensate density was small (< 10 13 cm~ 3 ) 
and, therefore, mean-field effects could be neglected. Af- 
ter that, the lattice was accelerated to a final velocity 
uvb (n — 1,2, ...). Each time the condensate was accel- 
erated across the edge of the Brillouin zone, according to 
L-Z theory a fraction r (see Eqns. ^ and ^ underwent 
tunneling into the first excited band. In Fig. |[ the av- 
erage velocity of the condensate after the acceleration is 
shown as a function of a and Uq along with theoretical 
predictions using the L-Z tunneling probability. If the 
final lattice velocity is v B , one finds a final mean velocity 
u m = (1 — r)vB, whereas a straightforward generalization 
of this formula yields 



= v B (l/r-l)[l-(l-r) n ] 



(T) 



for a final lattice velocity of uvb- In this frac- 
tion r of the condensate undergoes L-Z tunneling each 
time the Bragg resonance is crossed, with a remaining 
fraction 1 — r being accelerated further. Note that this 
result is only independent of the lattice depth (except 
through the tunneling fraction r) if the final lattice ve- 
locity is an integer multiple of nvs- As can be seen from 
Fig. |3|, agreement with theory is good. Again, it should 
be stressed here that owing to the small condensate densi- 
ties, these measurements do not differ qualitatively from 
those using cold but uncondensed atoms (111 . 
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FIG. 9. Calculated dependence of the parameter C on the 
trap frequency vtrap for two different lattice configurations. 
Mean-field effects are far more important for the larger lat- 
tice constant d — 1.56 /im, which results in a C larger by a 
factor of ~ 16 with respect to the counter-propagating config- 
uration (d = 0.39 fjxa). In calculating C, we assumed a typical 
condensate number N to t = 10 4 and used the Thomas-Fermi 
expression for the condensate peak density. 



VII. MEAN-FIELD EFFECTS 

A. Theoretical considerations 

In Bose-Einstein condensates, interactions between the 
constituent atoms are responsible for the non-linear be- 
haviour of the BEC and can lead to interesting phenom- 
ena such as solitons [|| and four-wave mixing with matter 



waves p0[ . As the atoms are extremely cold, collisions 
between them can be treated by considering only s-wave 
scattering, which is described by the scattering length 
a s . For rubidium the atomic mean-field interaction is 
repulsive corresponding to a positive scattering length 
a s = 5.4 ran. 

For a BEC in an optical lattice, one expects an ef- 
fect due to the mean-field interaction similar to the one 
responsible for determining the shape of a condensate 
in the Thomas-Fermi limit: The interplay between the 
confining potential and the density-dependent mean-field 
energy leads to a modified ground state that reflects the 
strength of the mean-field interaction. Applied to a BEC 
in a periodic potential, one expects the density modula- 
tion imposed on the condensate by the potential (higher 
density in potential troughs, lower density where the po- 
tential energy is high) to be modified in the presence of 
mean-field interactions. In particular, the tendency of 
the periodic potential to create a locally higher density 
where the potential energy of the lattice is low will be 
counteracted by the (repulsive) interaction energy that 
rises as the local density increases. 

The nonlinear interaction of the condensate inside an 
optical lattice with lattice constant d — ir/ sin(#/2)fc£ 
may be described through a dimensionless parameter |lq] 



C 



Trn a s 



- sin- 



■■(9/2) 



E B > 



(8) 



with g defined in Eqn. (A.2), corresponding to the ratio 
of the nonlinear interaction term and the Bloch energy. 
The parameter C contains the peak condensate density 
no (3lj, the scattering length, and the atomic mass to. 
In our notation, the parameter C always refers to the 
respective lattice geometries with angle 9. From the de- 
pendence of C on the lattice angle 9 it follows that a small 
angle 9 (meaning a large lattice constant d) will result in 
a large interaction term C. In our experiment, creating 
a lattice with 9 = 29 deg (i.e. d — 1.56 um) allowed us 
to realize a value of C larger by a factor of more than 10 
with respect to Rcf. |3^| using a comparable condensate 
density. Figure ^ shows our estimates for the nonlinear 
interaction parameter C realized by varying the magnetic 
trap frequency (and, thereby, the density no). 

In Ref. (l8) , the authors derived an analytical expres- 
sion in the perturbative limit (assuming Uq <C Eb) for 
the effect of the mean-field interaction on the ground 
state of the condensate in the lattice. Starting from the 
Gross-Pitaevskii equation for the condensate wavefunc- 
tion in a one-dimensional optical lattice (i.e. a one di- 
mensional Hamiltonian equivalent to that of Eq. 
they found that by substituting the potential dept 
with an effective potential 



Al), 



U 



U, 



eff 



Uo 



1 + 4C" 



(9) 



the effect of the mean-field interaction could be approxi- 
mately accounted for. This reduction of the effective po- 
tential agrees with the intuitive picture of the back-action 
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on the periodic potential of the density modulation of the 
condensate imposed on it by the lattice potential. For re- 
pulsive interactions, this results in the effective potential 
being lowered with respect to the actual optical potential 
created by the lattice beams. 
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FIG. 10. L-Z tunneling for different values of the interac- 
tion parameter C. The tunneling fraction r was measured by 
accelerating the condensate to db in a vertical lattice with 
d — 1.2 /j,m, and C was varied by changing the trap frequency 
and hence the condensate density no. The straight lines are 
best linear fits to the data. 



In order further to test the validity of Eqn. [nj, we used 
two different lattice angles and varied the condensate 
density by changing the trap frequency fl35| . The effec- 
tive potential in each case was inferred from the tunneling 
probability r for a fixed lattice depth. The results of these 
measurements are shown in Fig. [LlJ along with the the- 
oretical predictions. Clearly, the reduction of U e ff with 
respect to the non-interacting limit is much larger for the 
small lattice angle, as expected from theory. The general 
behaviour of U e f / as a function of C is well reproduced 
by our results. Replacing C by » 2.5C in the formula 
for U e ff/Uo leads to much better agreement with the ex- 
perimental data, which is consistent with the results of 
Fig. 0. 
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B. Experimental results 

In order to measure the effective potential U e ff, we 
assumed that the perturbative treatment described above 
can be extended to define an effective band gap AE e ff 
at the edge of the Brillouin zone which for a particular 
interaction parameter C can be written as 



o.oo 



0.05 



AE. 



U, 



eff 



eff 



U 



2(1 + 4C)' 



(10) 



where AE is the value of the band gap in the absence 
of interactions [j33| . One can then derive the interaction 
parameter C indirectly by determining the effective band 
gap AE e ff from the L-Z tunneling rate r, using Eqns.fjJ) 
and (|5|) with the band gap AE e ff. 

Figure [To] shows the results of measurements of the L-Z 
tunneling rate for two different values of C . In this exper- 
iment, in contrast to those described thus far, the lattice 
was adiabatically ramped up with the magnetic trap still 
on in order to maintain the condensate at a constant den- 
sity. After accelerating the condensate to vb 1 the mag- 
netic trap and the lattice were both switched off and the 
fraction r that had undergone tunneling (i.e. the frac- 
tional population that appeared in the zero momentum 
class) was measured after a time-of- flight. The effective 
potentials could be derived from the slopes of the linear 
fits in Fig. [l^ and were found to be markedly different. 
Comparing the experimentally determined values for C 
with those calculated on the basis of the experimental 
parameters, we found that the experimental values were 
larger by about a factor of 2.5. In fact, we expected the 
predictions of |L8| to be only approximately valid (see 
discussion belowTT 
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FIG. 11. Dependence of the effective potential U a g on the 
interaction parameter C. The different symbols indicate the 
two (horizontal) lattice geometries, with the empty symbols 
corresponding to the counter-propagating case and the full 
ones to the angle configuration (d — 1.56 /im). Solid line: 
theoretical prediction of Choi and Niu's expression for t/ e //. 
The parameters in these experiments were a = 23.4 ms -2 
and Uo = 2.2 E rec for the counter-propagating lattice and 
a = 3.23 ms~ 2 and Uo = 5.7 E rec for the angle geometry. 



In spite of the good qualitative agreement, we should 
like to point out that the model of Choi and Niu |l^] only 
approximately describes our experiment, as it assumes an 
infinitely extended condensate and neglects the radial de- 
grees of freedom of the condensate in the one-dimensional 
lattice. Especially the finite extent of the condensate in 
our experiment (in which only 6 — 30 lattice sites were 
occupied by the condensate, see the Appendix) should 
lead to non-negligible corrections. Also, the analysis of 
Ref. HH assumes a uniform condensate density across 
the entire lattice, whereas in our experiment there was a 
pronounced density envelope over the 6 — 30 lattice wells 
occupied by the condensate. In the above comparison 
with theory we calculated C using the peak condensate 
density. 



VIII. CONCLUSIONS AND OUTLOOK 

We have presented experimental results on the adia- 
batic loading and subsequent coherent acceleration of a 
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Bose-Einstein condensate in a ID optical lattice. In the 
adiabatic acceleration limit we have observed Bloch os- 
cillations of the condensate mean velocity in the lattice 
reference frame, whereas for larger accelerations and/or 
smaller lattice depths L-Z tunneling out of the lowest 
band occurred. The experimentally observed variation of 
the L-Z tunneling rate with the condensate density has 
been related to the mean-field interaction in the conden- 
sate leading to a reduced effective potential. Agreement 
with recent theoretical results is satisfactory. 
A natural extension of our work on mean-field effects will 
consist in checking theoretical predictions concerning in- 
stabilities at the edge of the Brillouin zone |3(| and the 
possibility of creating bright solitons by exploiting the 
nonlinearity of the Gross-Pitaevskii equation, which can 
compensate the negative group velocity dispersion at the 
band edge §371 ]. 
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APPENDIX: AN ARRAY OF COUPLED WELLS - 
RELEVANT PARAMETERS 

The possibility of studying the dynamics of a Bose- 
Einstein condensate spread out coherently over a large 
number of wells of a periodic potential, bearing a close 
resemblance to an array of coupled Josephson junctions, 
has inspired a host of theoretical papers in the past few 
years. On the experimental side, phase fluctuations [p"l|| , 
Josephson oscillations |l5| and the Mott insulator tran- 
sition |l2"| have been investigated, invoking concepts and 
notations inherited from the physics of Josephson junc- 
tions. In order to facilitate the comparison of our work 
with these studies, in this Appendix we report the values 
pertinent to our experiment for the various parameters 
that are important in the description of coherent quan- 
tum effects in an array of tunneling junctions. 

For the description of a condensate in an array of 
coupled potentials wells, the physical parameters needed 
to describe the dynamics of the system are the on-site 
interaction Ec, and the tunneling energy Ej. These 
quantities are defined in a variety of ways in the litera- 
ture |j38|-ft0[ . Our calculations are based on a variational 
ansatz of the total Hamiltonian 



H tot = H +g\y(r)\* 



-h 2 . 



2m 



U sm 2 (k L -r)+g\t>(r)\ 2 , 
(Al) 



with the interaction parameter g given by 

g = A-n-h 2 a s /m (A2) 
and the wavefunction ^S{r) given by 

^(r) = J2Mx~nd)^N^T)e ii> ^ t U(y,z). (A3) 

n 

Here, N n (t) is the number of atoms at site n, d n (t) is a 
site-dependent phase, and 4>(y, z) is the part of the wave- 
function perpendicular to the lattice direction. Basing 
the variational ansatz for ipo on a Gaussian of the form 
i/jq{x) — tT i /2 1 7r i/4 exp[— ^(x/a) 2 ) |42|, we obtain a mini- 



mum energy wavefunction of width a which, expressed 
in units of the width Oh — ^{Ua/ Er^)^ 1 ^ in the har- 
monic approximation of the potential wells, satisfies the 
condition 



exp 



(A4) 



This equation can be solved numerically to yield a juu. 
We now define the quantities Ec and Ej as follows |tlj : 



(A5) 
(A6) 



In the expression for Ec, the ID interaction parameter 
5id is defined as 



Ec = N s g 1B J dxip (x) 4 , 
Ej = - / dxipo(x)H tp (x - d). 



1 



.9id = g- 



TTOyOz 



(A7) 



where <J VtZ are the Gaussian widths in the y and z direc- 
tions of the radial wavefunction 



1 



1 



(z/a z 



(A8) 



The ID coupling strength <?id is equivalent to that de- 
rived by Olshanii in the case of a cigar shaped atomic 
trap [|43[ . N s — N tot /n occ is the mean number of atoms 
per lattice site, with N tot = J2 n N n the total number 
of condensate atoms and n occ the number of lattice sites 
occupied, as defined below. 
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FIG. 12. Dependence of the on-site interaction energy Ec, 
in (a), and of the Josephson frequency vc, in (b), on the 
non-linear parameter C for a lattice with d = 1.56 /im, 
Uo = 5.6 E rec and N to t = 10 4 . For these parameters, the 
tunneling energy Ej = h x 8 Hz. 
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FIG. 13. Variation of Ej and Ec with lattice depth (lattice 



constant d — 1.56 /im, Ntot = 
non-linear parameter C — 0.17. 



30- 

25- 

8 20- 
o 

c 15- 



10 ) for a fixed value of the 



ID- 
S' 



d = 0.39 urn 
d = 1.56 urn 



— i — 

20 



— I — 
30 



_40 

V trap ( HZ 



~l — 
50 



60 



I 

70 



FIG. 14. Number of wells occupied by the condensate as 
a function of the magnetic trap frequency for two different 
lattice constants. For the large lattice constant, the number 
of wells is small, so that finite-size effects are expected to be 
important. 



As the maximum lattice depth we could experimen- 
tally achieve in the counter-propagating configuration 
was ss 2 E rec , we have only calculated Ec and Ej numeri- 
cally for a lattice in the angle-geometry with d — 1.56 /im, 
as the present model only gives reasonable values for 
Uo/E rec > 4. Figure |l2|(a) shows the dependence of 
the on-site interaction energy Ec as a function of the 
non-linear interaction parameter C for a constant lat- 
tice depth Uq = 5.6E rec . The Josephson frequency 
v c = \/EjEc/h as a function of C is shown in Fig. [l2| 
(b) . In Fig. [13], both Ej and E c are plotted as a function 
of the lattice depth. 

Finally, we briefly discuss the variation with v tra p of 
the number n occ of lattice sites occupied by the conden- 
sate. In a rough approximation, this number is given by 
the diameter of the condensate as calculated from the 
Thomas-Fermi limit divided by the lattice constant d. 
Pedri et al. have used a more refined model [^6| to de- 
rive the expression 



n occ = 1 + - 1 



1-Kmvt 



rap 



15 



Ntota-i 



mnVtrap 



1/5 



(A9) 



Figure |l4] shows n occ as a function of v trav for two 
different lattice geometries. It is clear from this plot that 
in the angle-geometry, the number of wells occupied (< 
10) is small and hence we expect finite-size effects to be 
particularly important in this configuration. 
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